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Abstract 



We propose a Newtonian theory of semiclassical gravity based on the 
GRWm collapse theory, which we call GRWmN. The theory essentially 
amounts to modifying the Schroedinger-Newton (SN) equations of stan- 
^^ ' dard semiclassical gravity so that the SN wavefunction undergoes random 

f^ , collapses according to the GRW law. The primary motivation for this the- 

ory is to suppress 'cat state' solutions that arise from macroscopic super- 
positions, thereby ensuring that a semiclassical gravity field is macroscop- 
ically sourced by a single mass density occupying a single space-time vol- 
ume. We apply the theory to the mass scale at which Giulini and Grofiardt 
have numerically found the standard SN equations to predict gravitational 
wavepacket localization, and examine the possibility of experimentally ob- 
serving the predicted GRWmN collapse on top of wavepacket localization. 
5_2] [ We also consider the possibility of a general relativistic extension of GR- 

C^ . WmN, and argue why, in our view, semiclassical gravity theories such as 

GRWmN are worthwhile to develop. 

1 Motivation 

The problem of how to consistently couple a classical gravitational field to quan- 
tized matter was first addressed by M0ller and Rosenfeld, who proposed the 
modified Einstein equation (also called the "M0ller-Rosenfeld" equation) 
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where < T„„i >=< ip\Tnm\tp >■ Using standard quantum mechanics, this turns 
out to be the only way to incorporate a quantum description of the right hand 
side of (1) while keeping the left hand side a classical field [T]. Such a formulation 
of semiclassical gravity has well-known difficulties, however, among them being 
the existence of 'cat state' solutions. To illustrate these cat state solutions, it is 
convenient to first take the Newtonian limit of (1). 

Making the approximations g„„ = r;„„ + /i„™, |r»™|/TOO = |r""|/p « 1, 
and V << c, (1) reduces to the semiclassical Newton- Poisson equation 

VV(x, t) = -ATTGm\ijj{x, i)p, (2) 

with solution V{x,t) = —G J ™u_^/| (fix', and ip satisfying the nonlinear 
integro-differential Schroedinger equation, 

-ihdtipix^t) = V^^'(.T,i)+ ml/fa;, i)V'(a;,i). (3) 
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The N-body generalizations (ignoring the interaction potential term for simplic- 
ity) are as follows: 
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ihdttlj{xi...XN,t) = -Y^t; — V'^-ip{xi...XN,t) + y^miV{xi,t)tp{xi...XN,t), 

(5) 
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V{x.„t) = -G>] / ^'V^^- ,7 ^' dx[...dx'j,. (6) 
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The coupled equations defined by (2)-(3) or (4)-(5) are known as the Schroedinger- 
Newton (SN) equations. They describe a physical world in which the wavefunc- 
tion in configuration space drives the dynamical evolution of a mass density field 
(or a set of N mass density fields in the N-system case) in 3-space, the evolv- 
ing mass density field(s) sources an evolving classical gravitational potential 
in 3-space, and this gravitational potential couples back to the wavefunction, 
thereby altering the dynamical evolution of the mass density field (i.e. the 
so-called gravitational 'back-reaction'). 

The SN equations can be shown to admit cat state solutions as follows. 
Elaborating the example by Ford [l], suppose we have a quantum state tp = 
-k= [01 + (f>2], where each state in the superposition corresponds to a macroscopic 
mass distribution in a distinct location (e.g. a f 000 kg mass occupying a volume 
located on the left or right side of a room). Inserting ?/> into (2) gives 



VV = -A-kG 



-T7\<f>l\^ + ^\h\^ , (7) 
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or the prediction of a semiclassical gravitational field which is an average of the 
fields due to the two distributions separately (in this case, the gravitational field 
is the sum effect of two 500 kg masses on opposite sides of the room). However, 
we would expect that an actual measurement of the gravitational field should 
correspond to a single 1000 kg mass density source occupying a single loca- 
tion, but in different locations in different measurement trials. Unfortunately, 
such a measurement outcome is not predicted by anything in the SN equations. 
Moreover, Page and Geilker's torson balance pendulum experiment has already 
disconfirmed the gravitational field predicted by (6) [1]. 

It should be remarked that incorporating the effects of quantum decoher- 
ence does not get rid of these cat states (for essentially the same reason that 
decoherence doesn't solve the measurement problem); all decoherence can do is 
ensure that 4'i{q) ■ 4'2i<l) ~ (i.e. 0i and 02 have disjoint supports in config- 
uration space) for all q = {xi, ...,xn) so that there are no interference terms 
contributing to the r.h.s. of (6). It should also be emphasized that the gravi- 
tational self-localization effect discovered in numerical simulations by Salzman 
and Carlip [5j, as well as Giulini and Grossardt [3], does not solve the cat state 
problem either - all the self-localization effect potentially does is ensure that 
each state in the superposition will localize separate, 500 kg mass distributions 
around their respective locations in 3-space. 

It seems then that semiclassical gravity based on a formulation of QM in 
which there is just the wavefunction evolving by (3) is not empirically adequate. 
We propose that one seemingly natural way to solve the cat state problem (which 
is the cause of the empirical inadequacy) is to develop a semiclassical gravity 
theory based on the GRWm collapse theory [5]. Accordingly, we develop a 
nonrelativistic version of such a theory in the next section. 

2 GRWmN 

For a single system, the "GRWm-Newton" (GRWmN) equations are defined as 
follows. Starting from the GRWm mass density field 

m{x,t) =m\ip{x,t)\'^, (8) 

we can use this as a source in the Newton-Poisson equation 

VV(a;, t) = -AnGmix, t), (9) 

where V{x,t) = —G J 7^_^',| rf'^rr'. This gravitational 'self-potential' couples 
back to the wavefunction via (3), but now the latter undergoes discrete and in- 
stantaneous intermittent collapses according to the GRW collapse law. That is, 
the collapse time T occurs randomly with constant rate per system of NX = A = 
10~^^s~^, where the post-collapse wavefunction ipT+ = ^i'm't\/r4't is obtained 



from the pre-collapse wavefunction V't- = Hmty^Tipt through multiphcation by 
the Gaussian function 

4>T+ix) = ^gix - X)1/Vt-(x), (10) 

where 

9i^) = (^^^^'^ (11) 

is the 3-D Gaussian function of width a = IO^^jti, and 

C = C{X)= ( j (fxg{x - X)\tl)T^{x)A (12) 

is the normahzation factor. The coUapse center X is chosen randomly with 
probabihty density p{x) = C(a;)^, and the spacetime locations of the collapses 
are given by the ordered pair (X^Tk) . Between collapses, the wavefunction just 
evolves by (3). The generalization to an N-body system is as follows. The N- 
body self-potential and wavefunction now satisfy (4) and (5) respectively, while 
(8), (10), and (12) generalize to 

N 

m{x,t) =^m, dx'^...dx'j^5'^{x', - X)\ij{x'^, ...,x'j^,t)\^ , (13) 



and 



and 



ll)T+{xi, ...,Xn) = -7^9{xi - XY/'^^T-{xi, ...,xn), (14) 



C^C{X)=(j dx[...dx'r,g{x[~X)\ijT-{x'i.::x'^)A , (15) 

where i is chosen randomly from 1 , ..., N. 

The equations of GRWmN for a single system say the following - a wavefunc- 
tion in 3-space, which evolves by (3) and undergoes the random collapse process 
in (10), drives the dynamical evolution of a mass-density field in 3-space via (8). 
When the wavefunction collapses, it localizes the mass density field around a 
randomly chosen point in 3-space, with width 10~^m, and the probability of 
the randomly chosen point is largest where the mod-squared of the uncoUapsed 
wavefunction is largest, as indicated by (15). This evolving mass density field 
also sources a gravitational potential in 3-space via (9) , and this potential cou- 
ples back to the wavefunction via (3), which in turn alters the evolution of the 
mass density field via (8) again. 

For N-systems, the wavefunction lives in configuration space M.^^ , evolves 
by (5), and undergoes the collapse process in (14); this wavefunction drives the 
dynamical evolution of N mass density fields in 3-space via (13) so that when 



the wavefunction collapses, it randomly localizes the mass density fields around 
randomly chosen (non-overlapping) points in 3-space, each of width 10~^m, and 
with probability density given by (15). As before, each of these mass density 
fields acts as a source for a gravitational potential in 3-space that couples back 
to the N-system wavefunction via (5) - (6), which in turn alters the evolution 
of the mass density fields via (13) again. 

Like the SN equations, the Schroedinger equation for the GRWmN wavefunc- 
tion also admits cat states, but because the GRWmN wavefunction undergoes 
random collapses according to (10) or (14), which scales with the number of 
systems, those cat states are not macroscopically observable. (Also, the gravi- 
tational field produced by a cat state for a single elementary particle is presum- 
ably far too weak to be experimentally measured.) For example, for a massive 
object composed of Avogadro's number of systems, the collapse rate is ~ 10''i. 
So the individual mass fields composing the massive object will be localized 
around definite points in space frequently enough to give the appearance of 
a macroscopic mass distribution occupying a particular volume of space. Re- 
turning then to the example of a 1000 kg mass, it is clear that the number of 
systems needed in practice to compose such a mass distribution would imply an 
astronomically faster collapse rate; and when such collapses take place via (13) 
- (14), the result will be the appearance of a single 1000 kg mass localized on 
either the left or right side of the room (assuming the collapse center X for each 
system can take a binary outcome - either the left or right side of the room) 
with equal frequency. Correspondingly, the gravitational field measured with 
a classical test particle will look like it is due to only one mass density source 
at one location. In this way, the gravitational field predicted by GRWmN is 
consistent with that observed in the Page-Geilker experiment, in contrast to 
standard semiclassical gravity. 

It should also be noted that, in addition to the GRW collapse process, 
the branches of the GRWm wavefunction can undergo the gravitational self- 
localization effect observed in numerical simulations of the SN equations for a 
free Gaussian wavepacket. In particular, Giulini and Grofiardt found that for 
m = lO^^u and initial width of 0.5 microns, a Gaussian wavepacket will undergo 
self-localization, reach a minimum width of 0.4 microns in 30,000 seconds, and 
disperse again thereafter [J. As it has been suggested [3l |4j [6] that molecular 
interferometry experiments with macromolecule clusters may eventually reach 
this mass scale, it is natural to ask if GRW collapse might also be observable at 
this mass scale and perhaps happen 'on top of the self-localization effect. If we 
make the generous assumption that in GRWmN a mass of lO^^w corresponds 
to 10^° systems of 1 u, this gives an approximate collapse rate of 10"^-!-, or 
lO^s for each collapse. In other words, to have a chance of observing the GRW 
collapse effect, we would have to maintain the coherence time of the wavepacket 
for at least 33 times longer than the timescale for self-localization to reach the 
minimum width. It remains to be seen whether technological advancements in 
molecular interferometry that allow for maintaining coherence times of 30,000 
seconds will also allow for maintaining coherence times of lO^s or greater. Even 
so, we note that if self-localization is not observed at the mass scale predicted 



by GRWmN, this will be sufficient to falsify GRWmN as a semiclassical theory 
of gravity. And if self-localization is observed, it would be strong evidence for 
GRWmN or some dynamical collapse variant of GRWmN. 

3 Discussion 

To summarize, we outlined the cat states problem of standard semiclassical 
gravity and presented a modification of Newtonian semiclassical gravity called 
GRWmN to solve the problem. In this way, we have a proof of principle that one 
can have an ontological formulation of semiclassical gravity that's empirically 
adequate in its description of macroscopic mass superpositions. In addition, 
we indicated how GRWmN could be experimentally testable with molecular 
interferometry experiments in the near future. 

To the extent that GRWmN is based on a version of GRW theory with 
a 'primitive ontology!^ (GRWm), one could ask if other GRW theories with 
primitive ontologies could be extended to semiclassical gravity. Along with 
GRWm, perhaps the best known GRW theory with a primitive ontology is GRWf 
[5], where f stands for the 'flash' ontology, i.e. an ontology in which matter is 
represented by the space-time locations of the collapsed wavefunction, or the 
set F = {{Xi,Ti), ..., {Xk,Tk), ...}. Obviously the flashes can't on their own be 
used as a source for a semiclassical gravitational field, as they are just space-time 
locations. One could perhaps modify GRWf so that at each flash location a point 
mass is spontaneously produced, which then acts as a source for a semiclassical 
gravitational potentiao However, such a theory would entail violation of mass 
conservation every time a flash is produced, making it seem rather contrived, 
even though logically possible. A more straightforward and natural way we can 
see to retain the flashes in GRWf and also extend it to a Newtonian theory 
of semiclassical gravity is to add the presence of a mass density field in space- 
time which, upon collapse, is localized around the flashes. But this is essentially 
what happens in GRWm, to the extent that the GRW collapse localizes the mass 
density fleld around the space-time locations of the collapsed wavefunction. It 
would then seem artificial and contrived to continue to insist that the flashes 
constitute the primitive ontology in this approach. So the more natural GRW 
approach to semiclassical gravity seems to suggest the mass density ontology 
as the more natural one to describe the physical world, even for cases when 
gravitational self-interaction effects can be neglected. 

Concerning the collapse mechanism, although we based our formulation of 
semiclassical gravity on the discrete and instantaneous GRW collapse, it seems 
entirely possible to also use a continuous and non-instantaneous collapse mech- 
anism such as CSL (Continuous Spontaneous Localization) [7], together with 
a mass density field ontology. We used the GRW collapse mechanism sim- 
ply because it is already used in GRWm, a theory which is known to match the 



^Primitive ontology is defined in [5] as just the "variables describing the distribution of 
matter in 4-dimensional space-time". 

^This possibility was suggested to us by T. Norsen [personal communication] . 



predictions of standard nonrelativistic quantum mechanics for all current exper- 
iments while also having a primitive ontology that can be naturally extended to 
semiclassical gravity. 

Although we restricted our theory to Newtonian semiclassical gravity, it can 
certainly be extended to the M0ller-Rosenfeld equation. As we will show in a 
forthcoming paper, what needs to change for this extension is the collapse law 
and the mass density field definition - rather than using the nonrelativistic GRW 
collapse law, we would instead use the relativistic GRW collapse law developed 
by Tumulka in [5], and the relativistic stress-energy-momentum tensor field 
proposed by Bedingham et al. in [5]. However, such a theory will have an 
inconsistency that also plagues standard semiclassical gravity - the covariant 
divergence of the r.h.s. of (1) is nonzero upon wavefunction collapse, while the 
covariant divergence of the l.h.s. (i.e. the Bianchi identitiy) is always zero. 
As Tumulka has noted [personal communication], this inconsistency may mean 
that equation (1) with a wavefunction undergoing the GRW evolution does not 
possess any solutions. This is a question that needs further researcHf], in which 
case, a GRW version of semiclassical gravity seems to fare no better here than 
standard semiclassical gravity. 

Along with this inconsistency, such a theory would also inherit the stability 
problem of standard semiclassical gravity (i.e. the fact that equation (1) is a 
fourth-order system means that some solutions have runaway behavior), and the 
formally divergent expectation value of Tnm [2 ■ But these latter problems seem 
to be manageable - Ford [5] has noted for standard semiclassical gravity that 
there exist adequate renormalization procedures for < Tnm >, and there exist 
reasonable proposals for solving the stability problem by reformulating equation 
(1) to eliminate the runaway solutions or by regarding the semiclassical gravity 
theory as valid only for spacetimes that pass a certain stability criterion. 

In light of all these problems that seem to arise with semiclassical gravity, we 
would like to anticipate a potential critic who might ask why one should care 
about developing a physically consistent formulation of semiclassical gravity, 
especially since many physicists ultimately want a full quantum theory of grav- 
ity which presumably won't have the problems of semiclassical gravity theory. 
First, we point out that semiclassical gravity (in its standard formulation), in 
spite of its problems, has been used to derive several key results in theoretical 
astrophysics - in particular, the Hawking effect, eternal cosmic inflation, cosmo- 
logical pair production, the creation of naked black hole singularities, traversable 
wormholes, and warp drive spacetimes [2llll). 

Second, unlike the predictions of most quantum gravity theories, many of 
these semiclassical predictions may be empirically testable in the near future. As 
examples, it has been proposed that a condensed matter analogue of the Hawk- 
ing effect may be observed in experiments using superfluids with supersonic 
flow velocities [T^] • The density perturbations in the CMB spectrum predicted 
by eternal cosmic inflation may also soon be tested with the Planck Satellite's 
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that can be given a collapse interpretation and also satisfies < Tnm >',m= 0, but it is not 
clear to us how this prescription can be extended to a GRW theory. 



mapping of the CMB power spectrum [13] . It has also been shown by Visser et 
al. that the prediction of cosmological pair production could be tested in a BEC 
that simulates a quantum field evolving on an expanding spacetime |14| . And of 
course, we have the prediction discussed earlier of gravitational self-localization 
from the SN equations, which may be testable in molecular interferometry ex- 
periments in the near future. Given the potential testability of the various 
astrophysical predictions of standard semiclassical gravity, it also seems to us 
worthwhile to re-analyze them with a GRW approach to semiclassical gravity 
to see if the GRW approach gives differing predictions that could be observed 
in said tests. In this regard, it has recently come to our attention [Struyve, 
personal communication] that Landau et al. |15) have constructed and applied 
CSL variants of GRWmN to a scalar quantum field undergoing cosmic infla- 
tion, and shown that such models leads to different CMB power spectrums than 
the standard, scale-invariant Harrison-Zel'dovich one. They have also shown 
that such collapse models are currently consistent with data on the CMB from 
WMAP. It remains to be seen if such deviations are also consistent with data 
on the CMB from the more recent Planck satellite. 

Third, one of the major approaches to quantum gravity, namely, canonical 
quantum gravity, actually has equation (1) as a prediction of its semiclassical 
limit [T]. So the semiclassical limit of canonical quantum gravity presumably 
inherits all the same problems as semiclassical gravity theory taken as funda- 
mental. If so, this may also suggest that canonical quantum gravity would 
benefit from a GRW-type modification. 

Fourth, it has been suggested that if gravity is emergent, then it would be 
misguided to try and quantize gravity [TBI El [HI US] j for example, there exist 
the many variants of Sakharov's induced gravity approach, which have equation 
(1) as the emergent description of the coupling between quantum theory and 
gravity |16j . Like with canonical quantum gravity, then, the induced gravity 
approaches could also benefit from a GRW-type modification. 

In sum, we believe semiclassical gravity is a very worthwhile approach for 
trying to consistently incorporate quantum mechanics and gravity. And we 
believe dynamical collapse versions of semiclassical gravity, such as GRWmN and 
its possible relativistic extensions, provide a more promising route to a consistent 
formulation of semiclassical gravity than does standard quantum mechanics. 
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